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POISSON STATISTICS FOR MATRIX ENSEMBLES AT LARGE 

TEMPERATURE 

FLORENT BENAYCH-GEORGES AND SANDRINE PEGHE 


Abstract. In this article, we consider ,5-ensembles, i.e. collections of particles with 
random positions on the real line having joint distribution 


1 

Z^) 


|A(A)|^ 


4 2^i = 


‘dA, 


in the regime where /3 —>■ 0 as X —>■ oo. We briefly describe the global regime and then 
consider the local regime. In the case where A/3 stays bounded, we prove that the local 
eigenvalue statistics, in the vicinity of any real number, are asymptotically to those of 
a Poisson point process. In the case where A/3 —> oo, we prove a partial result in this 
direction. 


1. Introduction 

General /3-ensembles are collections of particles with random positions on the real line 
with joint distribution 

( 1 ) 

where dA denotes the Lebesgue measure on A(A) := ni<i<j<Ar(-^j “ U is a 

2 

potential with enough growth at inhnity (like the Gaussian potential Vg{x) = and 
Z]si{l3) is a normalizing constant. Their study is initially motivated by some considerations 
from physics: the probability distribution can be viewed as the equilibrium measure of 
a one dimensional Goulomb gas, but they actually appear to be connected to a broad 
spectrum of mathematics and physics, such as random matrices, number theory, lattice 
gas theory, quantum mechanics and Selberg-type integrals. In the case where /3 = 1, 2 or 
4, the probability measure of (1) is the joint distribution of the eigenvalues of a random 
N X N matrix M with density proportional to on the space of respectively real 

symmetric, complex Hermitian or quaternionic Hermitian matrices (see e.g. [6]). Besides, 
it was proved by Dumitriu and Edelman in [15] that when V is the Gaussian potential, for 
any /3 > 0, the probability measure of (1) is the joint distribution of the eigenvalues of the 
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random N x N tridiagonal matrix 


H = 


(9i X, \ 

X 2 92 X^ 

X3 ••• ••• 


V 


Xn 

Xn 9n ) 


(2) 


where the are some iV(0,1) variables and for all i, Xi = \/Yi, with Yi distribnted thanks 
to the r((i — l)f3/2) law, everything being independent. For general potential V, there is 
no random matrix representation. 

In the classical cases {/3 = 1,2,4), the combinatorial strncture and repnlsive interaction 
has been well-nnderstood for a long time via the theory of determinantal or Pfaffian pro¬ 
cesses (see e.g. [6] for references). The nnderstanding of these asymptotic spectral statistics 
to the fnll class of parameters f3 > 0 has recently mobilized a lot of research. For general 
/3, despite the lack of strnctnre, some enormons progress has been accomplished recently. 
For hxed f3, a few resnlts are now known. First, it is known from [8] that the empirical 
eigenvalne distribntion of the rescaled matrix converges weakly as N ^ 00 to a. 

probability measnre which is the semi-circle distribntion in the case of Ganssian potential. 
The local eigenvalne statistics in the large iV-limit are also qnite well nnderstood. In the 
Ganssian setting, at the edge of the spectrnm, Ramirez, Rider and Virag have shown in [21] 
that the eigenvalnes of — 2a//VJ) converge in distribntion to those of the so-called 

stochastic Airy operator. In the bnlk of the spectrnm, the limiting spectral statistics are 
asymptotically defined in terms of the Sine-/3 process, which is again defined as the solntion 
of a stochastic eqnation by Valko and Virag in [26] . In particular the authors show that the 
Sine-/d, which is translation invariant, has a geometric description in terms of the Brownian 
carousel, a deterministic function of the Brownian motion in the hyperbolic plane. Some 
advances on /9-ensembles have also been made by Sosoe and Wang [22, 23] and Bao and 
Su in [7]. 

The question of universality for these statistics has now become an important matter of 
interest: some enormous progress has recently been accomplished by Bourgade, Erdos and 
Yau in [11, 12, 13, 10]. Therein the authors consider general /9-ensembles (when the poten¬ 
tial V is and regular, or, in the hrst papers, convex and analytic). Assuming that the 
limiting spectral distribution (which depends on V) is supported on a single interval, they 
prove that the limiting eigenvalue statistics at the edge of the spectrum are given by the 
/9-Tracy-Widom distribution. The universality in the bulk of the spectrum is also proved. 
Another point of view to tackle /9-ensembles and in particular the quantitative aspect of 
the repulsion between eigenvalues has been developed in particular by Allez, Bouchaud 
and Guionnet in [3, 1]. They show in particular that when /9 < 2, /9-ensembles can be seen 
as an Y-dimensional process whose evolution is a mixing of that of N independent real 
Brownian motions and of that of a /9-Dyson Brownian motion. 
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The scope of this article is to understand the spectral behavior, at microscopic scale, of 
/9-ensembles in the case where ,9 —)■ 0 and N ^ oo (so that (3 depends on the dimension 
N). At macroscopic scale, such ensembles have been considered recently by [25] (see also 
the close model studied in [2]). Therein it is proved that when (iN —>■ c for some constant 
c > 0, the scaled empirical eigenvalue distribution of converges to the spectral measure 
of a deterministic Jacobi matrix, the density of which is explicit. When /9iV —)■ cx), the 
limiting empirical eigenvalue distribution of converges to the semi-circle distribution. Local 
eigenvalue statistics have not been considered yet. 

We here also consider the regime where /9 —)■ 0 and N ^ oo, but study the local eigen¬ 
value statistics. In [18], Killip and Stoiciu have considered the same question for circular 
,9-ensembles. More precisely they study CMV matrices (which are discrete one-dimensional 
Dirac-type operators) with random decaying coefficients. For rapidly decreasing coeffi¬ 
cients, the eigenvalues have rigid spacing while in the case of slow decrease, the eigenvalues 
are distributed according to a Poisson process. More precisely, they prove that local eigen¬ 
value statistics of /9-circular ensembles when ,9 —)■ 0 are in the large N limit those of a 
Poisson process. 

For real-symmetric ensembles, the same question has recently been considered from a for¬ 
mal point of view. Indeed, in [4, 5], Allez and Dumaz considered the ,9 —?• 0 limit of the 
Sine-/1 process and of the ,9-Tracy-Widom distribution. The ,9 —)■ 0 limit of the Sine-/1 
is also considered by Leble and Sefaty in [19]. The approach used by [19] is based on 
approximation theory while [4] use the diffusion representation of the Sine-/9 process to 
consider the limit ,9 —)■ 0. 

One would expect again to prove that when /I —)■ 0 simultaneously to —)■ oo, the eigen¬ 
values in the vicinity of a point u in the bulk of the spectrum exhibit Poisson statistics. 
In this text, we prove that this is true when Nf3 stays bounded as iV —)■ oo. In the case 
where ,9 —)■ 0 but Nf3 —)■ oo, we have a partial result which formally implies the Poisson 
statistics in the bulk, but does not allow to get a complete proof. 

In Figure 1, we compare this result with numerical simulations, giving a numerical evidence 
of the fact that the Poisson approximation works well (but gets less accurate as jS grows). 

Notation. For u = u{N) and v = v{N) some sequences, 

M -C n u/v — > 0 ; u ^ v u/v —)■ 1. 

N^OO N^OO 

Acknowledgment. The authors thank Alice Guionnet for her contribution to the proof 
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Bulk spacing of /3-ensembles for /3=Y, N = 1000 



Figure 1. Numerical comparison of the spacings of the A^’s with the 
spacings in a Poisson point process: histogram, in blue, of the debiazed 

spacings (A^+i — Xk)^J^ — for 0.05N < k < .95N, with (3 = for several 
values of c G (—1, 0), compared to the density, in red, of the exponential law. The 
histogram of the spacings is drawn thanks to 10^ independent realizations of the 
distribution of (4) for N = 10^, (we ordered the Afc’s in an increasing way so that 
Xk+i — Xk is actually the spacing between two consecutive particles). We see that 
the Poisson approximation works well, but gets less accurate as f3 grows (note 
however that all simulations of the last row were made with /3 almost equal to 1). 

2. Statement oe results 


2.1. Presentation of the model. For any q;,/1 > 0 and any iV > 1, we define 

A(A)^e“^ ^‘=1 '^^dA, 


Zjv(a,/1) : — 

with A(A) := ni<i<i<iv W - 


'AeM^ 


( 3 ) 


Let us now consider an exchangeable family (Ai,..., Aat) of random variables with joint 
law 

• • • ) ^Aat) := y-c-^A(A)^e 2 Ei=i 

Zjv(Q!, pj 

with Z]sf{a,/3) the normalization constant defined at (3). 


( 4 ) 
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2.2. Tridiagonal model and relation between a and /9. Let 

(9x ^2 \ 

^2 92 -^3 


E = 


X. 


\ 


( 5 ) 


Xat 

Xt<! Qn ) 

where the gi^ are some N{0, 1) variables and for all i, Xi = y/Yl, with distributed thanks 
to the r((z — l)/9/2) law. We known, by [15] or Section 4.5 of [6], that is the joint 
law of the eigenvalues of H. 

Note that TrP is centered and oETriL^ = ^ go that for 

a ~ 1 + iV/3/2, (6) 

the empirical eigenvalue distribution of H has asymptotic hrst moments 0 and 1. 


2.3. Global and local regime for bounded Nfi. The following proposition gives the 
limit of the empirical distribution of the Aj’s. The probability measure /i.y in question here 
has been studied in [1, 25]. 

Proposition 2.1 (Global regime for bounded iV/9). Suppose that N(3 —)■ 27 > 0 as 
N ^ 00 and a —> 7 + 1. 


a) 


Then under the law , 


the random probability measure 



i=l 


converges in probability to an even probability measure on M, depending only on 
7 , with moments defined at (26), satisfying m 2 = 1 . 


b) The measure is absolutely continuous with respect to the Lebesgue measure on M, 
with a density f^^ that is bounded on any compact set, and satisfies, for all x > 0, 

7+1 J.2 

/i..y(]R\[—a:, x]) < C^- - 

where is a constant depending only on 7 . 

c) p.y depends continuously on 7 > 0, is egual to N{0, 1) if and only if ■y = 0, and 
tends to the semicircle law with support [— 2 , 2 ] as 7 —t 00 . 

d) For each k, Xa.i{N-^ Y.f=i >^i) = 0{N-^). 
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The following theorem gives the limit local behavior of the Aj’s. 

Theorem 2.2 (Poisson limit for bounded N(3). Suppose that N(3 — > 27 > 0 as —>■ +cx) 
and that a ~ Nf3/2 + 1. Fix E eM.. As N ^ 00 , the point process 

N 

1 = 1 

with (Ai,..., Atv) ~ ! converges in distribution to the law of a Poisson point process 

with intensity 6 da; on M, for 

with is in Proposition 2.1. 

Remark 2.3. Note that the formula of 9 given at (7) should agree with the density of 
at E as given in [1, 25], but we were not able to prove it so far. 



2.4. Case where (3 N Using the fact that N‘^f3 3> N, one can easily adapt the 
proof of the following theorem from [6] . 


Theorem 2.4 (LDP hr (3 N ^). Suppose that as N ^ 00 , a = a{N) and {3 = I3{N) 
are such that a ~ N(3/2, N(3 3> 1 and (3 is bounded. Then for (Ai,..., Xn) distributed 
according to , the seguence of random probability measures Lj^ := A satisfies 

a LDP in the set of probability measure on M endowed with the weak topology with speed 
N^(3 and good rate function I defined by 


Hl) ■= jj f{x,y)p^‘^{dxdy) - 

with / : —)■ M U {+C)o} is the function defined by 

fix,y) = ^ - ^\og\x-y\. 

Moreover, the unigue minimum of I is achieved at the semicircle law 


a 


:= - a;2l|a.|<2da; 


27r 


and we have 


lim — 

N^oo m/3 


logZjt{a,l3) = - // f(x,y)cr^^-(dxdy) = 


( 8 ) 

(9) 


( 10 ) 


In the case where N(3 S> 1, as far as the local regime is concerned, we only have the 
following partial result, inspired from Johansson’s work in [17]. Below, we explain how 
formally, it allows to prove the convergence of local statistics to the ones of a Poisson point 
process and to identify its density. 







POISSON STATISTICS FOR MATRIX ENSEMBLES 


7 


Theorem 2.5. Let = P{N) and a = a{N) positive such that 

’ W-2a«l. (11) 

Let h : M —)■ M 6e a bounded function having 9 continuous hounded derivatives and 
(Ai,..., Atv) be distributed according to . Then as N ^ oo, we have 

log J ... ,dAAi) - A^/3 j h{t)da{t) —)■ J h{t)du{t) (12) 

for 

^ a ('^-2 + < 52 ) — l|x |<2 - / dx (13) 

Remark 2.6. The measure u of (13) is a classic correction to the semi-circle law (see e.g. 
[17, Rem. 2.5] or, more recently, [16]). 


Let us now explain how, on the formal level. Theorem 2.5 gives, for any E G 
the convergence of the point process 

N 

'^dNiXi-E) 

i=l 

to a Poisson point process with density 6dx on M, for 


(- 2 , 2 ), 


9 : = 




x^ 


271 


The hrst thing one has to notice is that 


9 = — exp / log|P-t|di/(t) 

277 J 

for u being as in (13) (the proof goes along the same lines as [ 8 , Lem. 2.7]). 
To prove it rigorously, we would need to prove that for 

R[^\xi,...,Xk) := J e^A:tiEf=i''iog|E+j^-A,|p^N-fc)^^^^^ _ _ _ 

for any Xi,... ,Xk G M, 

P) 

and that we have an upper bound of the type of (64). 

First, it can be proved (see Section 5.1) that as iV —>■ cx), for any hxed k, 

ZN-k{a,f^) ^ 

15) 


(14) 


(15) 


(16) 


(17) 


r\j 


27! 


(18) 
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Moreover, there is a universal positive constant M independent of N such that uniformly 
on N, k, 




ZN-k{0i, P) 

ZN{a,l3) 


< 



k 


(19) 


Theorem 2.5 can be rewritten as follows: for each hxed h as in the theorem, 




E (jv) = exp < iV/3 / h{t)da{t) + / h(f)di/(f) + £Ar(h) > , 


( 20 ) 


with £Af(h) 1. By (20), cutting on the right thanks to Lemma 5.1 and making as if 
the function : A h-)■ d- — X\ were C® (and close enough to the function 

h : A I—>■ fclog \E — A|), we should have 




for u as in (13). But by [8] p. 529, we know that for any E G (—2, 2), 
j k-^h{t)da{t) = j log \E - t\da{t) = ^ - 
so that we should have 

rI^\xi, ... ,Xk) ~ expi^NfS^^ -+ j h{t)du{t) + eN{h)^ . 

Besides, by (18), we have 

ZN-k{a,P) 

ZN{a,/3) \ 271 J ' 

Puting together (21), (22) and the fact that a ~ we should have 




exp < A^/3 / hAr(f)dcr(f) + / hAr(t)di/(t) + e:Ar(/i7v) 


exp < A^/3 / h{t)da{t) + / h{t)di'{t) + eN{h) 


( 21 ) 


( 22 ) 


ZN-k{a,(3) _ 
Zn{(^i Id) 


2 li? . . . ,Xk) 



k 

e 




^')+J h{t)du{t) 


= e'^ 


(23) 


with 6 as in (15). 
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3. Proof of Proposition 2.1 


Let H be as in (5). We shall prove that for any k, N tends to ruk, that the 


rukS satisfy Carleman’s criterion and that d) holds. By Skorohod’s representation theorem 
(see e.g. [14, Th. 2.3.2]) and a diagonal extraction, it will imply a). Part c) will be clear 
from the proof. Note hrst that if P is a r(t)-distributed variable, then for all k > 0, 



(24) 


Let us £x fc > 0, M G [0,1], let i = i{N) be such that i/N —> u and compute 
We have 



e 


where the sum is taken over paths £ : {0 ,..., fc} —?• {1, ..., n} such that 

• e:(0) = e{k) = i, 

• for all £ = 1 ,..., fc, e{i) — e{£ — 1) = —1, 0 or 1, in which case we say that i belongs 
respectively to D(£), F(£) or U(£). 

Note hrst that for such a path e, 


# D(£) = # U(£), # D(£) + # F(£) + # U(£) = k. 


(25) 


For any j G {1,..., A^}, we introduce 


F,(£) := = 

U,(t) := {f € U(£) ;£(«)= j} 


Then one can easily see, using (24), that 



j 


# Fj (e) is even ^ 


IS even 


I\ —^CXJ 



n ^ lui'ju + 1) ■ ■ ■ ( 7 M + #Uj(e:) - 1)) 


3 
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By the dominated convergence theorem, it follows that N ^ E Tr converges to 

+ [ JJ{#Fj(e:)!! X 7M(7M + l)---(7 M + #Uj(e:)-l)}dM, (26) 

£ Jn=0 j 

where the sum runs over paths £ : { 0 ,..., fc} —)■ Z whose steps are in {— 1 , 0 , 1 }, such that 
e( 0 ) = e{k) = 0 and for all k, jfFj{e) is even. 

Note that j^Fj{e) = #F(e:), whose parity is the one of k by (25), so that when k is 
odd, E(Fr^)jj = 0. Using (25) again, we see that when k is even, for any e, for any u, 

n ^ 7w(7w + 1) ■ ■ ■ (7W + #Uj(e:) - 1)) < max{l, 

j 

so that the m^’s satisfy Carleman’s criterion. It follows that the m^’s are the moments of 
a unique measure which depends continuously on 7 . Besides, d) follows from the fact 
that Cov((i7^)jj, = 0 as soon as \j — i\> 2k. 

If 7 = 0, then the only way for the term associated to £ in (26) to be non zero is that k is 
even and £ is the constant path equal to i. This proves that /xq = N^(0,1). The reciprocal 
is obvious, as the fact that /r.y tends to the semicircle law when 7 —)■ cx) (using the formula 
of the moments of the semicircle law in terms of Dyck paths, as in [ 6 ] or [20]). 

To prove the first part of b), we use Lemma 4.2 below. For any a G M and £ G (0,1], we 
have 

2 

/i((a — a + e)) < liminf Pa^g{\Xi — a| < e) < Cee^. 

V—>-00 

The second part of b) is a direct consequence of Lemma 4.3 below. 


4. Prooe oe Theorem 2.2 


4.1. Correlation functions. To prove the theorem, according to Proposition 6.1, we 
introduce the correlation functions of the point process ^N{\i-E), given by the formulas 


:= 7 .,.... 

ZN[a,p) [N - ky. jv 2 


with 


(27) 




First of all, we know that 


N-’^m 


< 1 


( 29 ) 
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and that as iV —)■ oo, for each hxed k, 

N-’^m ^ 

[N — k)\ N^oo 

Besides, for any M > 0, for any k > 1 and any xi,..., G [—M, M], we have 

A{xu...,Xkf ^ ({2MY\^ 


'-k<N- 


Bk(k-l) 

N 2 


\ N 2 / 

and as soon as /d (logiV)“^, for any fixed k and any fixed Xi,.. .Xk, 

A(xi,.. .,Xk)^ 




1 . 


(30) 


(31) 


(32) 


4.2. Partition functions. We know, by [6, Cor. 2.5.9], that 


Z,(a,n = 


Hence 


= aia4(^(^-i)-(N-fc)(A-fe-i)) -j-j- f27r)-^^2 ^ 2 ) 


ZN{a, /3) 


= aia4(2fc^-^(fc+b) 


-fc)(A-fe-l)) JJ" ^2 

j=N—k^l 
N 

n 

j=N—k-\-l 


r(i + f) 


+ i 


r(i + f) 


(33) 


By hypothesis, N(3 is bonnded and so is a. Let C > 1 be such that N[3 + a < C. Then 
we have, uniformly in k, 


ZN_k{a,/3) ^ / 0^*2^^ max[i,c] r \ 

Zi^{a,l3) ~ y^/^ minfi^ciT J 

Besides, as N(3 —> 27 > 0, for each hxed k, 

ZN-k{a,f3) , f (7 + 1 )"+^ y 

ZN{(y,(3) Yv^r(7 + i)y 


(34) 


(35) 


4.3. Uniform upper-bound on the correlation functions. 

Lemma 4.1. Let 1C be a compact subset of M. There is a constant C depending only 
on 1C and on the upper bounds on the sequences N(3, a such that for all /c, N and all 
Xi,..., Xfc G /C, we have 


lk<NR^^\xi ,... ,Xfc) < 
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Proof. Note that by (36) and (37), for each i G {1,... fc} and j G {1,..., iV — k}, we have 

n- 


^ ,8 ^ n-r^ . ^ r.a r+ N )^ + 


l-E + ^ - Xj\^ < {\E + ^1 + < 2^exp{/3 ^ 

Hence for C a constant (that might change from line to line) as in the statement of the 
lemma, 

W\xi,...,Xk) < 2^^^exp{^^(E+^)2} /e^^f=i"^'P^^“''^(dAi,...,dA7v_fc) 


2=1 


^ ^N-kjo^ jX) ^k^kCb'^ 


ZN-k{<X, /3) 

,kCb‘^ 

where we used (33) and the fact that for any x G [0,1/2], (1 — x)“^ < H 




□ 


Hence by (29), (31), (34) and the previous lemma, we have proved that b) of Proposition 
6.1 is satisfied. It remains to prove a) for 6 given by (7). 


4.4. Preliminary estimates. 

Lemma 4.2 (Bulk eigenvalues). There is a constant C depending only on the upper bounds 
on the sequences a and N(3 such that for any a G M and e G (0,1], 




Na§_ 2 


Proof. We have 

y{N) 


P^y{\X,-a\<e) = -- 


- / dAie-.'^^ 

P) 2Aig[a±£] J 


./ueKJv-i 


Note that if |Ai — a| < e, the for any j, |Ai — /ij| < £ + |/Xj — a\. Moreover, for all x G M, 
f > 0 , we have |x| < , hence 


|x| < 2 e 16 . 

{x + y)'^ < 2 x^ + 2y‘^, 


Using also the fact that 
we get that if |Ai — a| < £, 

|Ai - /ijl^ < (^ + \Pj - a\)^ < P 


(36) 

(37) 


Hence 


2 ^^e 




<y(|Ai-a|<5) < 

Z7v(a;, 


I dAie 2N / A(/i)^e h 4^ 


2 opfjpp 


M/i 
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We have 


Hence 


,2 


/i,- - /3- 


Aa — 13, (3a . 2 

[lij + 




Aa — (3 Aa — (3 


PS(\\,-a\<e) < 


Aa Aa 

-——-- / dAie 2'^l / A{fi)^e ^ 4a iMj + 4a-/3-> 

ZN{a,(3) JAiSlaie] i/xSR^-l 

= ^N-i(a - (3/A, (3) 

ZN{a,^) 

Na0 ^2 

< Cee^o—/3)“ , 

we we used (33) to upper bound partition functions quotient. □ 

Lemma 4.3 (Largest eigenvalues). There is a constant C depending only on the upper 
bounds on the sequences a and N(3 such that for all x > 0, 


pTAM > 0 <c 


g-(40-/3) 2|- 


X 


(38) 


Proof. We have 

p5’(|AiI>:^) = z 

ZN[a,P) J|Ai|>a 

Note that by (36) and (37), 

|Ai — fij\^ < (|Ai| + l/ijl)^ < 2^e^ 8 


— f dAie 2^1 [ A{pfYl\^i - 2Ea^"d/i 

^P) J\Xi\>x , 


Hence 


2N0 


dAie 




Zj\j-(^a, (3') 

2^/5 Ziv-i(a-/?/4,/3) 


a (3/A Z]if (^a(3^) J\\l\>x^ya—|3/4 

Then we conclude using (33) and the fact that for all y > 0, 

P + OO 


A(/i)^e 2 (“ 
e 2 


_H , e 2 

e 2 dt < - 

y 


(39) 

□ 

Lemma 4.4 (Tail of the empirical spectral law). There are some constants C, c depending 
only on the upper bounds on the sequences N(3, a such that for all N and all x > 0, 

+. 3- X% 


P. 




N 


> x) < Ce~ 
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Proof. We use again the tridiagonal matrix model of (5) for of the Aj’s. We know that 
+ ■ ■ ■ + A^ has the same law as Trif^, with H the matrix introduced as (5). Note that 


a 


by the well known convolution relations between Gamma-distributed variables, — Tr 
has a r((^)-distribution for := A^(l -\- {N — l)f3)/2. Hence 


P 


+.< P(G>5f^). 


a./3 jy 

Then, one concludes using the concentration inequalities for Gamma variables (see [9] p. 
28-29) which say that for all u > 0, 

P{G>ip{l + u)) < 

□ 

Lemma 4.5. For C, c as in the previous lemma, for any u, M, 6 such that cM"^ > 26 > 0, 


1 < / < 1-f 


a,l3 


2C6eP^ 

cM‘^-29' 


Proof. The integral above rewrites 


i + e 


|Ai - u\ 
M 


V 1) > x)dx. 


Now, note that as log(|A — m| V 1) < (A — u)^. 


pW(^ 


1 


5^1og( 


|Ai - u\ 
M 


VI) > U < 


(Ai — -|- 


-F {Xn - uf 


(40) 


N 


> M'^x) 


— a,l3 \ 

^ p(N)Al + .+ A 

- W./3 ^ N 

Then one concludes using (40) and the previous lemma. 


< + 2^^ +.+ 2A^ -F 2u‘^ ^ ^2^^ 


N 


N 


> 


M^x 


-u\ 


□ 


4.5. Convergence of the correlation functions. Let us now prove a) of Proposition 
6.1 for 9 given by (7). Note hrst that by b) of Proposition 2.1, we know that 


log |P — a;|d/i..),(a:) < oo. 


(41) 


Besides, by (27), (30), (32) and (35), it suffices to prove that for each k and each 
xi,... ,Xk G M, the quantity R^^\xi, ..., Xk) defined at (28) satisfies, as N ^ oo, 


rI^\xi, ... ,Xk) —t exp {2'jk I log\E - x\dfi-y{x) 

N^oo ' 


( 42 ) 
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4.5.1. Upper-bound. Let us prove that for any fixed k and Xi,...,xt, 

..... 

N^oo t J 

For e: > 0, set 

(<"' < isr-'’ 

suffices to prove that for any e: > 0 small enough, we have 


(43) 


As R-l < RI and, by (41), inf / log(|x — E\ y e)dfr-f{x) = / log |x — E\dfr-f{x), it 


hmsupi?^^’^^(xi,... ,Xfc) < exp^2'yk j log(|i? — x| V e)d/x.y(x) |> . (44) 

N^OO 


Note now that for any M > 0 large enough, as 

I E ~\- X'! TV — X 

log(|F^ + Xi/N — x| V e) = log{(|i? + Xi/N — x| V e) A M} + log ( ^ V 1 


M 


the function rI^’^\xi, ..., Xk) rewrites 


Q0 Eti EfLi'' iog{(|i?+^-A,|V£)AA/}g^Eti Ef.i" iog('''+^ LI, 


M 




4 ELi EfEi" log{(|i?+^-A, |V£)AM} p{N-k) 


./3ELi EfEi" log{(|i?+^-A,|V£)AM} 


F;;7"'(dAi,...,dA^_fc) + 

/ Eti EfEi" iog(^^±%Fh VI 


l|Ay''’(dA„...,dA, 


:=X 


■“v"” 

:=Y 


By Proposition 2.1, we know that under the law the L°°-bounded sequence of 

k N-k 

random variables fd EE log{(|P + ^ —AjI Ve) AM} converges in probability, as —>■ oo, 

i=l j=l 

to 


2'yk J log{(|P — x| V e) A M}d/i.y(x) 

(one gets rid of the ^’s by noticing, for example, that the convolution of probability 
measures is continuous with respect to the weak topology and that Sx./n converges to ^o)- 


Note that by choosing M large enough, one can make J log{(|P — x| V e) A M}dfi-y{x) 
as close as we want from Jlog(|P — x| V e:)d/i.y(x). Moreover, one can easily adapt the 
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proof of Lemma 4.1 to see that f ^^(dAi,..., dA 7 v-fc) is bounded by a constant 

independent of M, hence by Cauchy-Schwartz, to prove (44), it suffices to prove that 
f F^P^^~^^(dAi,..., dAiv-fc) can be made as small as we want if M is large enough. Note 
that for 


log 




Yi := e 

we have Y = Yi ■ ■ ■ — 1. Thus by the Holder inequality, it is enough to prove that for 

M large enough, earh Y) can have its fc-th and 2fc-th moment as close as we want from 1, 
which is a direct consequence of Lemma 4.5. 


4.5.2. Lower hound. To obtain the analogous lower bound 

> exp| 27 A; [ log |P - x|d/i.^(x) , 

N^oo 

we observe that hrst by Jensen’s inequality and then by exchangeability, 

k N-k 

,{N-k)i 


(45) 


s ^ log|B+^ - A,|P),)“'(dAi,...,dAK_i) 

= {N-k)P /'j]log|B+^-A,|Fiy‘'(dA„....dAM_ 7 .) 

1=1 


Hence as {N — k)(3 —> 27 and the triplet {N, (3{N),a{N)) satishes the same hypotheses 
as {N — k,(3{N), a{N)), it suffices to prove that for any hxed x, we have 


X 


liminf / log|P+^ - Ai|Pf^^(dAi,...,dAAr) > / log |P - t|d/r-^(f). 


N^OO 


N 


As, by exchangeability, /i..,, is also the weak limit of the distribution of Ai under , we 
know that for any £ > 0, 


X 


liminf / log(|P +- All Ve)P]^^(dAi,... ,dAAr) > / log |P - t|d/i.^(f) 


N^OO 


N 


(and one can get rid of ^ for the same reason as in Section 4.5.1 above). Hence it suffices 
that for e small enough. 


lim sup 

N^OO 


(logdE + £ - All V £) - log IE + - All) Pi5(dA„... .dAw) 


can be made as small as desired. But for any random variable X > 0, 

E[log(X V £) - log(X)] = E[(log£ - log(X))lx<£] = [ dt. 


Here, by Lemma 4.2, there is a constant C such that 

X 

N 


PS{\E+^-Xi\<t) < CL 
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which allows to get the desired bound. 


5. Case where N~^ -C /S < l/log(A^) 


5.1. Partition functions: proofs of (18) and (19). It follows from (33) that for all N, 

for all 1 < fc < iV, 



1 


r(l + (iV-£)/3/2) 


e=o 



r(i + (iv - £)/?/ 2 ) 


But by the Stirling formula, for z —>■ +oo, 


r(l + z) ^ \/271 z exp {z log( 2 ;) — z} . 


(46) 


Note that our hypothesis on (3 implies that for any hxed i = 0,... ,k — 1, we have 


{N - i)^/2 log((iV - £)/3/2) -{N- e)/3/2 = N^/2 \og{Np/2) - N^/2 + o(l), 


so that, as (11) implies that N^ log ^ <C 1 and (3 log a <C 1, 




By (11), we get (18). The upper bound (19) comes in the same way, noticing that the 
error in (46) is uniformly bounded on ^ > 0. 

5.2. Tail estimate. 

Lemma 5.1. Let f3 = /3{N) and a = a{N) be satisfying (11). Then there is a constant C 
depending only on the seguences a and (3 such that for all N, for all x > 0, 



(47) 


Proof. We have 




|Ai|>a: .//xeK^-l 


dAie 2'^! f A(/i)^ |Ai —/ijl^e 


Note that by (36) and (37), 


|Ai - < (|Ai| + \ij.j\f < 2^e^ 's ' 
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Hence 


P 5 ’(|Ai|>a:) < 


2^/3 


Zn(,0',(3) i|Ai|>x 
2^/3 


dAie 




I ^leR^- 




J d/i 


. -- / _e 2 ^ dAi X Z7v-i(a — y, /?) 

^A(a,/3)\/a! - /?/4 y|Ai|>*y'a-/3/4 4 

x(a —/d/4) ZN{a,f3) 

x{a-l3/i) ' Zj,(a,^) 


< 


Let us now use (19) and for example the fact that {1 — y) ^ < e^^ when y G [0,1/2], We 
get 

f5’(|A.I>a^) < ^^exp{-^(a-P/i) + N^py(Sa) + NP/(2a) + Nl3/2}, 
which allows to conclude, as we already noticed that (11) implies that 2a ~ iV/3 S> 1. □ 


5.3. Proof of Theorem 2.5. We hrst dehne the probability measure on 

N N 


pN 


■= ^h7^7^^(^)^exp{-^^x2 + /?^h(a;i)}da:i---dx/v, 




i=l 


2=1 




^a,y 


where 


Z';^{a,(3) = j 

is the normalisation constant. Let, for z = 1, 2, 


(48) 


(49) 


/ p^^^’(xi,...,x;v)dx,+i---dxAf 

be the z-th correlation function of 

Lemma 5.2. Let ip : M C be a function on M such that the real and imaginary parts 
of tp' are bounded below. Then we have 

1^ ^ 2 JJ + J {I3h'{t) - at)^p{t)u]p'^{t)dt 

+ f ^p'{t)u]p’^{t)dt = 0. (50) 

Jr 
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Proof. As (50) is linear in -0, one can suppose 0 to be real-valued. Then for 6 > 0 small 
enough, the function y + 6ijj{y) is an homeomorphism on M, hence one can make the change 
of variable Xi = yi + 0f){yi) in (49). We get 

= [ n Ivj - Vi + ( 51 ) 

JyeR^ l<j<j<Ar 

N 

e-f JJ(i + 0^\y.))dy^ ■■■dyN 

i=l 


Let us compute the derivative, with respect to 9, ai 9 = 0, of the RHT of (51). We have 

de,e=o\yj -yi + Oi'ipiyj) - f){yi))f = - yif, 

yj yi 

we also have 

N 


2=1 


and 


Hence 




N 

P'^h'{yi)ij{yi), 

2=1 


N N 

00,e=o JJ(1+ ^0'(i/i)) = ^0'(2/i). 
2=1 2=1 


06 ,6=0 (RHT of (51)) 


^ ^ ^ • ■ ■, yN)dyi ■■■dyN 

i<i<j<N yj yi 

N 

/ yii^{yi)p%{yi,---,yN)dyi---dyN 

^ f 

/ h\yi)'ip{yi)p%{yi,...,yN)dyi---dyN 

N 

+ y] / 'ip\yi)p%(,yi,---,yN)dyi---dyN 


We get exactly (50). 


□ 
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Now, we define, for 2 : such that $ 52 : > 0, 
Hn(z) := [ 


Um(z) ; = 


1 l,h 


and 


Kn{z) := N 


_ z-t ^ ^ J z-t 

kN{s,t) := Nu]f'{s)u]:l'{t) - {N - 

1 1 1 


(t)dt 


{z — t){z — s) 2{z — ty 2{z — s)^ 

1 ff (t - sf 

2JJ {z-t)‘^{z-sy 


kN{s,t)dtds (52) 


-kp^{s, t)dsdt. 


We also introduce 


U{z) := -{z - Vz'^ - 4) 


(^ 5 ( 2 ;) > 0), 


(53) 


where when 3fJ(;2) > 0 (resp. < 0) , y/z"^ — 4 is computed with the determination of the 
square root on C\(—cx),0) (resp. on C\(0,+cx))) with positive values on [0,+cx)) (resp. 
such that y/—l = i). It is well known that U is the Stieltjes transform of the semicircle 
law a. 

Lemma 5.3. On the upper half-plane, we have 

N/3{Un-U){2U-^z + Um-U) = Am := /3N-^Km- 2/3+ {2-/3)U'^ + {2a-N/3)U\ 

(54) 

Proof. We shall apply the previous lemma with ■^(t) = -. Note that _ 


z — t 


t — s 


{z-t){z-s) 


, so that we have : 


N{N-l) II ^^^K f{s,t)dsdt 


t — s 


N{N — 1 ) m ^^( s , t)dsdt 

{—Nk]sf{s,t) + N‘^u]/l^{s)u]/f^{t))dsdt 

/cAr(s, t)dsdt 

1 1 


{z-t){z-s) 
1 

{z-t){z-s) 
= N^UN{zf-N 

= N^UM{zf-N 

-N 


{z-t){z-s) 
1 


{z — t){z — s) 2{z — tY 2{z — sf 


kj^is, t)dsdt 


- --^A;Ar(s,t)dsdf 

{z t) 

N^UMizf-KM{z)+nU'M{z) 
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(where we use the fact that for any function/(t), jj f{t)k]s[{s,t)dsdt= / f{t)u]^^{t)dt). 
We also have 


t'ip{t)u]^^{t)dt = 1 - zUn{z) ; / = -U'j^{z), 




so 2N ^(3 ^ X (50) rewrites 

U^{z)^ - N-^K^(z) + N-V„{z) + 2N-'H^(z) + ^(1 - zUn(z)) - -^U'„{z) = 0 


z.e. 


+ ^ = - 2iv-‘ir«(.-) + (^ - ^)C2;vW. (55) 

One gets (54), using the well known equation f/(x)^—xl7(x) + l = 0 (see [6, Eq. (2.4.6)]). □ 


A key step in the proof of the theorem will be to prove that as iV —)■ cx), 

(3Kn{z) < N. 


(56) 


We shall now prove (56). Let {yi, - ■ ■ ,yN) be a random vector with distribution 
and for g G C;,(M, C), dehne the random variable 


N 

hjv(^) := '^9{yi), 

i=l 


with variance T,]\f{g) := E,[\jl%{g) — ¥.[jl’]^{g)]\‘^]. As 

K^{z) =E[{fi%{g) -E[i,%{gm 


for g{t) = [z — t) \ we have 

\Kn{z)\<En{^). 

Z L 


(57) 

(58) 


Note that as 

m%{9)\^] 


N 

E[g{yi)g{yj)\ 

Li=i 

A^(A^-l) JJ g{s)g{t)u^j^'^{s,t)dsdt + N 


u 


l,h 

N 
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we deduce 

^N{g) = N{N-1) jj -u]:l'{s)u]f{t))dsdt 

+N (^J u^^^{t)\g{t)\^dt - \ j u]:l'{t)g{t)dt\‘^ 

Lemma 5.4. There is L > d and c > 0 such that for any fixed function g, 
Sjv(5') < N{N-1)[[ g{s)g{t){u%^{s,t) - u];f^{s)u]^^{t))dsdt 


-L,L]2 

L pL 

u\\p,f+\\‘^A+ I / odT 


+N{ I ^uT{tMt)\^dt 


uf,^{t)g{t)dtf +4N 


’-L 


7-2||„|| g-cY/3 


Proof. Using the fact that | — -^ll^lloo, we see that the probability measure 

dehned at (48) and its normalization constant can be controlled thanks to the 

probability measure and its normalization constant: for any Borel set A C we 
have 

It follows that up to a change of the constant C, Lemma 5.1 is also true for P^p^\ which 


allows to conclude. 


□ 


This lemma allows to reduce the problem to a compact set, and after rescaling, one can 
turn the compact set in question to [—1/2,1/2] : we deduce, as in [17], that for w := z/L 
and 

p%{t, s) := u%^{Lt, Ls), p]^{t) := u^A^{Lt). 

we have 

Tn{ -) < N{N - 1) / - Gn{x ^ - - ){t)dt + 4N^e-^^^{^z)-fi (59) 


' z — t 


J-l/2 W — t ' W — X' 

where Gjq is the operator on L^([—1/2,1/2], dx) defined by 
GnUW) = I f{s){p%{t,s) - p]^{t)p]^{s))ds + f y/W - f pN{s)f{s)ds | . 


iV-1 \ L‘ 


'- 1/2 


Thus to prove the estimate of interest (56), we have to upper bound: 


(60) 


Following [17], we introduce the integral operator on L^([—1/2,1/2], dx) with kernel 
Pw{t, s) = • Then (s, t) \—)■ :^^Gn{x ^)(s) is an integral operator satisfying 
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the hypothesis of Theorem 2.12 of [24]. This trace class operator is nothing but 
thus by this theorem, we have 

r" —GN{=^mdt = Tt{P^Gm). (61) 

W-S 

We will not here recall all the arguments used in [17] Lemma 3.12 and Proposition 3.9 to 
estimate this trace. The proof transfers to our setting using minor modihcations (essentially 
replacing N by N/3). Note that the important feature of h is that is Lipschitz on compact 
sets in our case. Thus we simply state the hnal estimate we will use in this article, namely 
the following lemma ; 

Lemma 5.5. We have TT{Pu,G]sr) < GN~^ log(iV) for some constant G. 


It follows, by (58), (59) and (61), that 

\PKj^{z)\ < /3S^(^—) < GNp\og{N) + +4Ar2/5e-'=^^(3z)-2. 
z — t 

This is of course N, so the estimate of interest (56) is proved. 


As /31og(A^) <C 1 and N(3 — 2a -C 1, by (54), we deduce that, uniformly on compact 
subsets of C"*", one has that 


NPiU^iz) - U{z)) 


2U'{z) 
2U{z) - z 


z _ 1 

^ 2-4 -4 


(62) 


One recognizes easily that the RHT of (62) is the Stieltjes transform of the null mass 
signed measure z/ of (13) (to do that, use the fact that U{z), given by (53), is the Stieltjes 
transform of the semi-circle law and then use an integration by parts). 


The rest of the proof of the theorem is an easy adaptation of the proof of Theorem 2.4 
in [17] (p. 169-172). The main idea is to dehne 

P(A) ;= logEe^^^f-i'^^^^) - N(3X j h{t)da{f), 

to notice that 

dxF{\) = N [ h{ty/\t)dt, 

Jr 

to prove (12) for derivatives and to deduce (12) by dominated convergence. We use (62), 
namely the convergence 

in Fourier transform manipulations, precisely via the formula 


SN(t)dt (Tyz > 0) 
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with (5jv(t) := 


6. Appendix: Poisson limit eor point processes 

Let A be a locally compact Polish space and /x be a Radon measure on A. We consider 
an exchangeable random vector (Ai ,... ,Xn) taking values on A implicitly depending on 
N, with density with respect to /i®^. We dehne, for 1 < k < N, the k-th correlation 
function on A^ by the formula 

:= [ p^^\xi,...,XN)dp‘^^-^{xk+i,...,XN). 

Proposition 6.1. Suppose that there is 9 > 0 independent of N such that the correlation 
functions satisfy: 

a) For each k > 1, on A^, we have the pointwise convergence 

Ri^\x,,...,Xk) (63) 

N^oo 


b) For each compact /C C A, there is Q/c such that for all k, N, on IC^, we have 

ik<NRi''\xi,...,Xk) < e’k (64) 

Then the point process converges in distribution to a Poisson point process with 

intensity 6dp as N ^ oo. 


Proof. Note that the Poisson point process M with intensity 6dp is characterized, among 
random random Radon measures on A, by the fact that for any compactly supported 
continuous function / on A, we have 




= exp 


9 


(qRR — l)dp{x) 
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So let us fix / a compactly supported continuous function on X. Then, with the convention 

= 1 . 


N 


= Ejj(l + (e' 


,/(Ad 


1 )) 


2=1 


PC{l,...,n} ieP 


N k 


k=0 'a a 

N k 

k=0 i=l 


This proves the proposition. 


□ 
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